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|I described possible alternative combustion catalyst configurations and 
detailed analysis of relevant performance factors, i.e., mass-transfer limited hydrocarbon 
conversions, pressure drop, catalyst temperatures, post combustion conversion of hydrocarbons 
and CO, etc. This report represents a comprehensive collection of my findings and comments. 
Attached are the final MathSojft Mathcad 6 analyses of the suggested alternative combustion catalyst 
configurations. 

Alternative Combustion Catalyst Configurations 

Two approaches are suggested for the catalytic combustion of gasoline vapors in air. 
These approaches represent more conventional and perhaps less expensive alternatives to the PCI 
Microlith™ technology. They are summarized as follows: 

• Monolith with bypass air A standard 3.5", corrugated FeCrAl alloy, 600 
cells/in 2 monolith should be able to reach « 86% (naphtha) to >99.8% (natural 
gas) fuel conversions. Using an air by-pass would allow longer residence time 
and higher adiabatic temperature rise (which is probably necessary for natural 
gas but may not be for diesel fuel or naphtha) within the combustor to ensure 
fast catalytic reaction rates. Higher temperatures within and at the exit of the 
catalyst monolith would also increase homogeneous combustion reactions, and 
help convert residual fuel vapors in the region immediately downstream of the 
catalyst but upstream of the turbine. A controller for bypass flow could help 
with turn-down, startup, and tuning steady-state conversions. 

• Radial flow catalyst bed Provided that the flow rate can be reduced to 
< 3 m/s, a thin packed bed of catalyst beads or granules may also achieve 
conversions similar to those of the monolith and PCI Microlith™ catalyst screens 
at equivalent pressure drop. The bed might intrude upon your fuel injection and 
premixing zone, but this could actually benefit combustion (within limits, e.g., 
±30% fiiel/air ratio) by providing locally higher adiabatic temperatures. A 
packed bed has fairly good thermal conductivity and is virtually immune to 
thermal shock problems. A screened frame and possibly metal plates (or 
alternatively for robustness, silica-coated silicon carbide) could be located within 
the bed to guide the flow. Downstream hydrocarbon and CO conversion also 
would be necessary for low emissions as in the case of the monolith. 



• Exhaust converter - A 400 cells/in 2 post combustion monolith located just 
ahead of the regenerator exhaust stream inlet could provide a very significant 
reduction in hydrocarbon emissions regardless of the configuration and 
performance of the combustion catalyst. An emission control monolith should 
remove about 99% of the residual unburned fuel components and > 99.9% of 
the CO. 
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Analysis of Catalyst Performance 

Mathsoft's Mathcad 6.0 Plus Professional program was used to perform detailed 
calculations of pressure drop and conversion of C\, C 5 , C 7 , and C9 hydrocarbons (methane, 
isopentane, cycloheptane and nonene) and CO. The analyses generally assume catalytic rates are 
fast under Capstone's specified combustion conditions. The analyses use various catalyst 
parameters such as particle/channel size and bed length for the two upstream (monolith and fixed 
bed) configurations. Only monoliths were considered for the downstream 1-atm conversion of 
residual CO and hydrocarbons. 

Gas Properties 

Initially the Mathcad document describes and defines various gas component and mixture 
properties and functions (pp. 1-13) that are used in the subsequent calculations. The gas and vapor 
components are identified with their physical parameters (p. 1) and the composition of combustion 
mixtures as a function of fuel conversion are defined (pp. 2-3). Mixture densities, diffusivities of 
individual fuel component in the mixtures, viscosity, and thermal conductivity are then defined and 
plotted as functions of composition and temperature (pp. 4-9). The composition of various fuel 
mixtures is shown to minimally effect mixture physical properties (< 0.5% under such lean 
conditions). These properties are then redefined as empirical functions of temperature for 
simplification of subsequent analysis (pp. 6-9). Heat capacities, heats of combustion and enthalpy 
functions of the mixtures are similarly defined (pp. 9-11), and the initial composition necessary for 



an adiabatic temperature rise of 51 1°F are determined for each of the four fuel components (p. 1 1 
and re-entered on p. 2). Finally, transport functions, e.g., Prandtl Number, Schmidt Number, 
thermal diffusivity, momentum diffusivity, etc., are defined (pp. 12-13) prior to definitions of the 
catalyst configurations, geometry, and flow conditions. 

Catalyst Configurations 

The physical parameters of several catalyst configurations (honeycomb monoliths and 
packed beds) and the appropriate flow conditions are described (p. 14) and used in the subsequent 
analyses. 

Combustion Monolith Performance 

The pressure drop values for (standard) 3.5 or 4.0 inch length 400 and 600 cells/in 2 
monoliths are determined as functions of flow rate (pp. 15-16, and see figure on p. 16). Heat and 
mass Nusselt functions are defined and used to calculate single-pass fuel conversion vs. monolith 
O depth (p. 17-19). It is possible to determine laminar flow hydrothermal Nusselt functions for 

l JB square or sinusoidal channels (which may be more representative of actual monolith channel 

shapes), but this was not included in the current analysis. Accordingly the Nusselt functions 
hi cannot be considered more accurate than about ± 30%. 

Q1 ' The results of the performance analysis for the monoliths are summarized in a table (p. 18) 

j~] for the two cell densities, the two lengths, and 0, 20, and 40% by-pass air. The transport-limiting 

conversion of fuel components is plotted as a function of pressure drop for constant flow rate but 
Q extended monolith length (p. 19). High fuel conversions (= 99.5%) are possible with the 

80 600 cells/in 2 monolith at 1% pressure drop (but it requires 10 inches length). 

ru 

fU Post Combustion Converter 

P Analysis of post-combustion downstream ( 1 -atm vs. 3-atm upstream) conversion of CO 

Is- 

™ and fuels and the pressure drop in 400 and 600 cpi monoliths shows excellent results for such a 

wide area monolith (pp. 20-25). After defining component diffusivity functions and describing 
monolith configuration (pp. 20-21), the pressure drop and conversions of CO and fuel vapors are 
calculated and plotted as functions of monolith depth (pp. 21-23). Finally, the results of an 
analysis of the overall conversions of individual hydrocarbons and CO including upstream 
combustor, homogeneous combustion, and the exhaust converter are tabulated (p. 23). 

Differential (Local) Nusselt Analysis 

A different analysis, where the local hydrothermal Nusselt functions are defined and 
integrated over the length of a monolith (pp. 24-25), are used to check the integral formulas. 



These values (p. 25), which are derived from different sources (Kays and Crawford, "Convection 
Heat and Mass Transfer," 3rd, 1993, vs. Rosner, 'Transport Processes in Chemically Reacting 
Systems", 1986) show higher conversions than the integrated values, e.g., 94% conversion for 
C 7 H U vs. 86% with the integral expression. This difference appears to be caused by the greater 
effect of hydrothermal entrance (as opposed to hydrodynamically developed flow with thermal 
entrance) and shows the integral formulas to be conservative by a factor of about two. 



Combustion with Packed Beds 

Pressure drop and fuel conversion functions are defined and used to analyze the 
performance of packed beds (pp. 26-3 1). Pressure drop is a strong function of the superficial fl ow 
rate, nearly second order (pressure drops are calculated and plotted against superficial linear flow 
velocity for various sizes of catalyst particles on p. 27). Gas to particle Nusselt numbers for 
q packed beds are developed plotted for methane, isopentane, cycloheptane and nonene for three 

5 sizes of catalyst granules (pp. 27-28). Shape factors are not included in this analysis, so the 

S results represent values for packed spheres with equivalent effect on AP, not with equivalent 

diameter. The conversion of fuel components (again assumng very fast catalytic rates) is 
pi determined (pp. 28-29) for the oblique flow geometry (which is defined on p. 16). 

0^ For full perpendicular flow such as with an annular cylindrical packed bed rather than the 

^ conical bed , the pressure drop function is much greater for equivalent conversion of fuel 

q components. This is caused by the higher power dependence of AP on flow velocities (form 

m drag), while the conversion per unit depth of bed increases almost linearly with increasing flow 

[U rate. At flows higher than about 2 m/s the catalyst granules begin to act as bluff bodies which 

2 cause the flow to accelerate around them. While this increases backmixing which helps 

U convection, it causes a greater effect on the pressure drop. Conversions for methane and 

cycloheptane are plotted as a function of pressure drop for beds of i-mm and 1/8-in. granules (pp. 
29-31). For oblique flow (2.5 m/s superficial velocity through the conical bed) with the small 
granules, C 7 H 14 conversions are calculated to be very high (> 99.99%) at 1% pressure drop, 
while for perpendicular flow (8.2 m/s) with the 1/8-in. pellets, C 7 H 14 conversions are calculated 
to be very low (~ 62%) at 2% pressure drop. Clearly some radial configuration must be used to 
achieve high conversions with the packed beds without a pressure drop penalty. 

Required Catalytic Rates 

Finally, the specific catalytic rate constants required for conversions of methane and 
cycloheptane in a 600 cells/in 2 monolith are determined (pp. 32-35) for Capstone's conditions. 
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The effect of gas diffusional restrictions within the macroporous (pores between agglomerated 
washcoat particles within a catalyst layer) and microporous (within the porous particles) is defined 
for CH 4 and C 7 H 14 (pp. 32-33). The effect of rate constant on surface temperature (for a point at 
the entry to the channel) is determined (p. 34) and shows the rates necessary to move fully into the 
transport limited regime, i.e., when the concentration of fuel near the surface approaches zero and 
temperature approaches a limit. This analysis could be turned around by calculating the conversion 
for a given catalytic rate constant. 

But, further analysis of the expected conversions requires measurements of the actual rate 
of the washcoated catalysts and catalyst beds. The analysis concludes with a complex rate 
expression for methane combustion by supported PdO that was based on measured rates (in the 
absence of transport limitations). This rate was based on PdO, an hexa-aluminate support 
(LaAl, j0 18 ) with a much lower surface area than the La-stabilized supported used in the Capstone 
monoliths and pellets. The interception of the curve showing the surface temperature as a function 
of catalytic rate constant (as on p. 35) and the curve showing the actual rate constant of the catalyst 
as a function of catalyst temperature (plotted on p. 37) given the actual performance of the catalyst. 




More detailed calculations would integrate the conversion through the monolith or packed 
bed and it would include the thermal conduction along the monolith channel or within the bed. 
This would represent a oneKlimensional calculation (axial surface and gas phase temperatures) 
which is only an approximation to the actual two- or three-dimensional convective transport and 
reaction problem. Such calculations could provide estimates of extinction conditions for specific 
catalysts. Such calculations and further analyses are beyond the scope of the current comments, 
but could be performed by extending the current Mathcad 6 program. 



Conclusions 



Catalyst Performance 

Several conclusions can be drawn about the transport limiting performance of a catalysts 
used to combustion gasoline and methane fuels for Capstone's 24 kW microturbine: 

I ) Commercially available small cell (600 cells/in 2 ) honeycomb monoliths 
should be seriously considered as an alternative to the PCI Microlith 
technology. A catalytic monolith configured as annular flow 600 cells/in 2 
honeycomb can achieve very high conversion of methane (99.99%) and 



modest conversion for gasoline (94% for cycloheptane) etc. with low 
pressure drop. 

2) If the fuel air mixing can be moved upstream sufficiently to allow use of a 
longer catalyst monoliths, significantly lower emission could result without 
unacceptable (> 1%) pressure drop ratio. For 1% AP/P drop (with greater 
monolith length), transport limited conversion of C7H14 should reach 99.5% 
in an active 600 cells/in 2 honeycomb monolith. 

3) Post combustion monolith catalysts should be considered for use in EZEV 
turbines, regardless of the method of fuel conversion during combustion. A 
short 400 cells/in 2 Alpha-IV metal monolith coated with a layer of supported 
Pt located at the inlet to the recuperator can give very significant emission 
reductions for CO and all hydrocarbons because of the lower linear velocities, 
high exhaust temperature, greater diffusivities at 1-atia Conversions of CO 
and CH4 can be very high, while conversions with C7H14 and other gasoline 
fuels is less, it still can bring the gasoline range fuel components in the 
exhaust to a few ppm with the assumption that only a modest degree of 
homogeneous combustion occurs upstream of the turbine. 

4) Provided that significant homogeneous combustion occurs within the 
channels and downstream of the monolith (« 80% conversion giving 60 ppm 
C7Hi4 and 840 ppm CO leaving the turbine) and that a 400 cells/in^post 
combustion hydrocarbon oxidation catalyst was used downstream, gasoline 
emissions should fall below 0.5 ppm. 

5) The oblique flow/conical catalyst packed bed could also be considered as an 
option to the Microlith™ and honeycomb monolith catalysts. The 
performance of the packed beds can be very good in terms of fuel conversion 
with the smaller catalyst granules. Yet, these beds must have some element of 
radial flow as the pressure drop for a flow linear flow rate will cause 
excessive pressure drop. A word of caution here, packed beds are more 
susceptible to flash back caused by the higher degree of backmixing and 
regions where local residence times can be quite long relative to the screens 
and monoliths. 



Catalyst Durability, Robustness, and Costs 

^^^^^^^^^^^^^^fcconcems were expressed about the tolerance of 
comBusuoi^ataly^ high fuel/air ratios and local 

unmixedness. While clearly high degrees of F/A control and uniformly 
premixed combustion gases are desired, the combustion catalysts should 
show robustness sufficient to survive short excursions of off-specification 
F/A ratio conditions as well as showing the durability necessary for long 
active life under normal operation at high temperatures. 

2) The greater the thermal mass of the catalyst, the greater its robustness in the 
face of locally or globally high F/A. In this regard the packed beds would be 



superior to the fine screens of the Microlith™ and the thin walled metal 
monoliths. 

3) A packed bed also has greater ability to disperse locally high F/A mixtures to a 
greater extent than the Microlith™ and to a much greater than the honeycomb 
channels. High thermal conductivity also improves robustness in locally high 
F/A (hot spot) conditions. Here the packed bed is perhaps superior to the 
metal monolith because of it greater radial conductivity, while the screens 
have very little ability (radiation) to disperse heat from locally hot zones. 

4) The Alpha-IV alloy honeycomb should be much more resistant to fast 
overheating to high temperatures than the PCI microlith. The Alzeta/Grace 
electrically heated monolith previously examined by Capstone (Nomac) had 
entrance channels too large to give high conversion except with long length 
and smaller channels downstream, hence the large volume and weight. 

5) Packed beds represent the least cost catalyst for combustion because they use 
configurations and materials common with commercial VOC incineration 
catalysts. The metal monoliths have low costs because of their widespread 
acceptance and increasing use in automobile catalytic converters and as VOC 
catalysts. The premixed, gasoline-fueled, and recuperated nature of the 
combustion process for Capstone's microturbines is a good fit (with a slight 
push to higher pressure, higher throughput, low pressure drop, and greater 
robustness) with the commercial VOC catalysts available from several 
sources, including PrototechAJnited Catalysts. It would not appear necessary 
to develop new catalyst and catalytic engineered combustors for Capstone* 
microturbine combustion applications. 



The following table summarizes the findings of this analysis: 



Property 




Ranking 




Fuel Conversion 


Packed Beds > 


Screens > 


Monoliths 


Pressure Drop 


Monoliths > 


Screens > 


Packed Beds 


Durability 


Packed Beds ~ 


Monoliths ~ 


Screens 


Robustness 


Packed Beds > 


Monoliths > 


Screens 


Cost & Availability 


Packed Beds > 


Monoliths > 


Screens 



Final Recommendations 



Use of Catalysts in Existing Engines 

The following recommendations are made on the basis of improving the existing 
microturbine catalytic combustion performance: 

1) TnstalMQO cells/in? l-atm monolith with a Pt-based catalyst for post combustion 
CO and hydxocaiboiu^Qyal^This unit should provide the greatest improvement 
in emissions wiffTfitSe effectoh pressure drop at modest cost (» $200/unit 
installed). It could be used with all catalytic and even with non-catalytic burners 
to reduce CO and hydrocarbon emissions. However, it cannot be effective in 
reducing NOx emissions from non-catalytic burners. 

2) Use a 600 cells/in 2 Alpha-IV alloy metal monolith of 4-in. length with premixed 
fuel and air as a lower pressure drop alternative to the PCI Microlith™ system. 
Fuel conversion may be lower than the PCI unit, but with lower pressure drop; 
at equivalent AP/P, the monolith should show superior fuel conversion. A small 
amount of by-pass air could help with turn-down and could be used to improve 
emissions during normal operation by increasing residence time and allowing 
higher temperatures in the homogeneous zones within the channels and 
immediately downstream of the combustor. 

Use of Catalysts in Advanced Engines 

For advanced combustion systems, the fuel air mixing could be moved upstream of its 
current position to allow great premixing and to give greater volume for the catalyst. 

1) Nominally 1-mra to 1/16-in. packed beds in a conical configuration with about 
2-3 m/s linear superficial bed velocity could give good performance, especially 
with a small bypass air stream for post catalyst hydrocarbon reduction. An inconel 
screen could be located upstream to suppress flashback and ceramic foams or 
honeycomb monoliths could be used as mechanical supports to increase robustness. 

2) A longer 600 cells/in 2 Alpha-IV alloy metal monolith (from 3.5- to 8- or 10-in.) 
could also decrease hydrocarbon and CO emissions with acceptable AP/P by 
relocating the F/A mixing zone upstream. 



Additional Catalytic R&D Activities 

Additional Analysis Fuel combustion rates can be measured for specific catalysts 
under Capstone's combustion conditions. The rates can be incorporated in an extended analysis of 
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combustor performance for use in predicting fuel conversion, specifying required catalyst activity 
maintenance (e.g., Is it close to extinction at 1000°F inlet temperature?), and assisting in the design 
of current and next generation combustors (e.g., How much by-pass air is necessary? What are the 
limited on F/A ratios?, etc.). 

i 

Materials ^flHH|||[^||HHH^IIIilHiifllBi^HHHH 

MHMpnore robust units could be developed with ceramic structural components and state-of-the-art 
catalysts capable of resisting short (a few seconds) temperature excursions to 3000°F (1650°C). 
The attrition resistance and crush strength of the Norton S A 6576 catalyst pellets supplied to 
Capstone by SRI International also need improvement (their very large specific volume of 
macropores, while desirable, leads to poor mechanical strength), and other catalysts (e.g., 1.5-mm 
extruded, Pt-promoted, Mn-substituted hexa-aluminates) could be examined for activity and 
attrition resistance. 

the longer-term.^HIIHIHmHiHHili^lHli^l^^lHHik^ 
HHlfcSRI could coat steels with chromium and aluminum (in situ preparation of FeCrAl alloy 
films) to prepare inexpensive, high-performance recuperators and monoliths for use in advanced 
higher temperature engines. 

Systems - Additional engineering analysis and catalyst preparation and testing could 
provide Capstone with cost effective assistance in the design of EZEV turbogenerators with 
alternative fuels such as LPG, CNG, and diesel for hybrid vehicles. LPG could be used in the 
current designs with combustion and exhaust converter monoliths; CNG could be used with 
combustion and exhaust converter monoliths, but with by-pass air to increase F/A ratios and 
catalyst temperatures; and diesel fuel could require upstream mixing with long combustion 
monoliths or conical beds and longer exhaust converter monoliths. 





Gas Properties and Constants 



*****•***********+******#**** 



cm-0.01m J^kg m 2 sec 2 W^J-sec" 1 Pref - 1 atm 

N=kgmsec" 2 Pa=Nm' 2 atm = 101300 Pa Mm = 0.000001m 



J 



R =8.3 143 Rc = 8.3 143* Jmole 1 KT 1 

K 

n Gas MolWt a(U) c(LJ) 

1 CH4 MW Q = 16.044 gm a Q : = 3.785- 10" 10 m £ 0 ^148.6 K 

2 02 MW t = 31.999 gm o { =3.464 10" 10 m t { := 106.7 K 

3 N2 MW 2 =28.018-gm c 2 =3.798* 10" 10 *m ^ =71.4^ 
4C02 MW 3 =44.010 gm a 3 = 3.941* 10" 10 m := 195.2 K 



5H20 ,n 

MW 4 := 18.015 gm o 4 : = 2.641- 10" 10 m e 4 = 809.1 K 



Tb5 =301K Tc5 =460.4 K Pc5 =3.3910°Pa Vc5 =306 



h 6C5H12 Tc5 
UJ 2nnelhyW)ulwi0 e, :=_ ^ = V c5 3 10 10 m 



3 



p MW 5 =72.151-gm e 5 =365.59994*K o $ =6.73866* 10 10 'm 

ob 

[H Tc7 =572.2K Tb7 =374.1 K Pc7 := 3.47 1 0 6 Pa Vc7 =368 

rU 

r*j 7C7H14 j 
£ methyl^yclohexane g = _Tc7_ a . = Vc7 3. 10 -io. m 

6 1.2593 6 

MW fi =98.189 gm c 6 =7.1661' 10~ i0 'm e 6 =454.37942* K 

Tc9 := 592 K Tb9 = 420K Vc9 =580 Pc9 := 2.34- 10 6 Pa 

l 

8C9H18 Tc9 \ m 

? c, a, := Vc9 3 10 -m 

1-nonene ? t 2593 i 

MW ? := 126.243 gm a ? =8.33955- 10~ 10 m ^ =470. 10244 *K 

9 CO c g =91.7K o 8 =3.6910" 10 m 

MW g = 28.01 gm 
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Define constants for transport integrals: 



ajia 0 =1.16145 



0]Hl x =0.52487 



opb Q =0.14874 



apbj =0.7732 



aDa Q = 1.06036 



cDa, =0.193 



oDb Q =0.1561 



aDbj =0.47635 



=2.16178 onb 2 =2.43787 cDjl, = 1.03587 cDb 2 = 1.52996 



Define transport integrals: 



oDaj = 1.76474 



Qp(T):=c^ 0 -exp[-(am> 0 -lii(T))]+ J] <Vb^T) 

i^=l 



aDb 3 =3.89411 



fin(3) = 1.0394 



QD(T) =aDa 0 .exp[-(aDb 0 lii(T))]^ J] cDa^e^-oDb^l) 

1 



OD(3) =0.95002 
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Set Inlet Conditions: 

ych40 =0.01183 
y0c7h 14 =0.0021859 



P =3atm 

y0c5hl2 = 0.0028968 
y0c9hl8 = 0.0016878 



These values for yOi have been determined via enthalpy 
functions on p. 1 1 for temperature rise from 1 000F to 1 51 1 F. 



Define gas composition w/r methane using moist (2%) air: 



y(yo,x) 



yO(l-x) 
(l-y0)0.205-2ycti40x 
(l-y0)0.775 

(l-y0) 0.02+2yO*x 
0 
0 
0 



y(ych40,0.1) 



0.0106 

0.2002 

0.7658 

0.0012 

0.0221 

0 

0 

0 



Zy(ych40,0.5) = I 



CAPSTONE fv">00 



( Define function for calculating stoichiometric gas composition for combustion of 
2-methyf-butane using moist (2%) air: 



y5(y0,x) := 



(l-yO) 0.205- 8yOx 
(l-y0)-0.775 
5(y0-x) 
(l-yO)0.02 + 6yOx 
yO(l-x) 
0 
0 



1 + 2*y0-x 
Xy5(yOc51il2,0.5) = l 



y5(y0c5h!2,0.5) = 



0 

0.192262 

0.770523 

0.007221 

0.02855 

0.001444 

0 

0 
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Define function for calculating stoichiometric gas composition 
for combustion of methyl-cyctohexane using moist (2%) air: 



y y7(y0,x) := 



(1- y0) 0.205- 10.5 yOx 
(1-yO) 0.775 
7(y0x) 
(l-y0)0.02 + 7y0x 
0 

yO(l-x) 
0 



1 + 2.5yOx 



y7(y0c7hl4,0.5) = 



pj Zy7(y0c7hl4,0.5) = 1 

3 

H : Define function for calculating stoichiometric gas composition 
for combustion of 1-nonene using moist (2%) air: 



0 

0.19255 

0.771199 

0.00763 

0.027532 

0 

0.00109 
0 



y9(y0,x) 



(l-y0)0.205- 13.5 yOx 
( I - y0) 0.775 
9(y0x) 
(1 -y0)0.02 + 9y0x 
0 



yQ(i-x) 



1 + 3.5 yO x 



y9(y0c9hl8,0.5) = 



Zy9(y0c9hl8,0.5) =1 



*********•****•-** 



0 

0.192692 

0.771413 

0.007573 

0.02748 

0 

0 

8.414148- 10~ 



Define functions for evaluation of mixture properties (ideal gas) given fractional fuel conversion and temperature. 
Density and average molecular weight: 



pmix(x,T) '•- 



I 7 \ 

V yCydrtO.x^MW^ |~; MWmix(x) := 
\ kk = 0 



2 tfycMO.x^-MW, 
\kk = 0 



kk 



pmix(l,300*K) = 3. 470573* kg* m 



MWmix(l) =28.485 *gm 



Define diffusion coefficients for specific binary diffusion pairs: 



Difl(n,m,T) := 1.8829 10 

□ 

u 



i-27. 



1 \ 



MW 



+ MW ] 

m/ 



0.5 



o -i-a 

n nil Qj) 



•10 gm -cm sec K 



( £ n ^m) 



0.5 



P 

Pref 



/'gj 

^0^(0,2, 



300-K) =0.07417'sec 1 W 



Diff(2,7,300*K) =0.0l597 sec 1 *cm 2 



W Define functions for use with various fuel air mixtures: 

a 

l-y(ych40,x) 0 



" Dch4(x,T) 



ru 

□ 



Dc5hl2(x,T) = 



y, y(ych40,x) jgB 

. Difl(0,jgas,T) 
Jgas= 1 

1 - y5(y0c5hl2,x) 0 



Dc7hl4(x,T) = 



Dch4(0.5,298K) =7.31327- 10 6 -m 2 -sec 1 



Dc5hl2(0.5,298K) =2.28466*10 6 -m 2 sec 1 



' ySCyOcSh^.x).^ 

. 2_i Difl(5,jg«s,T) 
jgas= 1 

I - y7(y0c7hl4,x) 0 

y y7(yQc7hi4,x) jgas 

Difl(6,jgas,T) 

jgas= 1 

l-y9(y0c9hl8,x). 
Dc9hl8(x,T) = - 

y y9(y0c9hi8,x) jgM 

. £-> Difl(7,jgas,T) 
Jgas= 1 

****«« *********** ************************ A************************ *»#*♦** #***♦#*#*****♦ 



Dc7hl4(0.5,298 K) = 1.94254' 10 6 m 2 *sec 1 



Dc9hl8(0.5,298 K) = 1.53123M0 6 -m 2 'sec 1 



Display gas mixture properties 



idsp =0..9 



Tdsp idap =(750+ 50idsp)K 



Tin 



(1000- 32) - + 273.15 
9 



K Tin=810.92778 K 



Tout 



(1511- 32) + 273. 15 | K Tout = 1094.81667 K 



Plots (below) of Di for CH4, C5H12, etc. show large gap between CH4 and the other hydrocarbons 
and show very little effect of fuel conversion (even for methane) on Di because the fuels are so dilute. 

Therefore we can choose methyl cydohexane as the molecule representative of naphtha and ignore 
the effect of conversion on gas properties 
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o 



8*10 



6M0 * - 



%G — 



Vc5hl2(.5,T6ap. 



.4*10 J ~ 



idsp) 
idap) 
dip) 



no 




700 



800 



900 



1000 



1100 



1200 



Tds pid 



7.510 2 K 



8.5 lO^K 



910 2 K 



9.5- K^K 



1 lO^K 



1.05- lO 3 ^ 



LMO^K 



1.15 10 3 K 



l^lO^K 



3.61119 



10" 5 m 2 sec 



4.02644 10 5 m 2 sec 



4.45884 



10 5 m 2 sec 



4.90801 10" 5 m 2 sec 



5.37357 



10" 5 m 2 sec" 



5.85522 10" 5 m 2 sec" 



6.35265- 10" 5 m 2 sec' 



6.8656 



10' 5 m 2 sec" 



739384 



10" 5 m 2 sec" 



7.93712 



10" 5 -m 2 sec" 



1.001 



1.0005 



0.9995 - 



0.999 




700 



800 



900 



1000 



1100 



1200 



Tdm 



Tds pid» 


Dc7hl4(.5,Tdsp.J 


7.5- lO^K 




1.03475* 10* 5 -m 2 *sec l 


sio^k 




1.15717- 10" W-sec" 1 


8.5 lO^K 




1.28482 10" 5 -m 2 sec 1 


9 K^K 




1.41758- 10* 5 -m 2 -sec' 1 


9.5- lO^K 




1.55533 10" 5 m 2 sec l 


1 lO^K 




1.69798 10" 5 m 2 sec 1 


1.05- lO'-K 




1.84541- 10" 5 -m 2 -sec _1 


1.1 10 3 K 




1.99755- 10" 5 m 2 sec" 1 


1.15 lO^K 




2.15429 10" 5 m 2 sec" 1 


1.2 10*K 




2.31555 10" 5 m 2 sec" ! 



idsp 



^ Define interpolation formulas for Och4 and Oc7h14 vs. T(K) at x=0.5 via Jandel Scientific's Tablecurve program 

Dch4T(T) :=exp|- 11.718527 + 1.6257167 h||j - 4677 ^ 09 K j.^ Dch4(.5,100O K) =0.585522-sec"' W 

Dch4T(1000K) =0.585523'sec -1 'cm 2 



Dc7hl4T(T) =exp(- 12.923061 + 1.62605581n||j - 82 48 0p 98 K j .ggL Dc7hl4(.5, 1000K) = 0.169798-sec~' 'cm 

Dc7hl4T(1000-K)=0.169798sec~ l -cm 2 



Define viscosity functions for major constituents: 



□ 



2.6693 10 



26 



(MW 0 -T) 0J 1Q5 



cm 



iff 



I 



(mw 2 t) 0 - 5 

1.116(MW 4 T) 



2.6693- 10 



2.6693- 10 



-26 



10 5 



cm | / gm 

IT 



26 



10 3 . 



sec 



cm / gm 



2.6693- 10" 76 - 



10 s 



cm) /gm 
K 



0.5 



0.5 



0.5 



(3.11510" 8 cm) -Qui 



\514-k/ 



2.6693- 10" 



26 



10 s 



sec 



cm \ j gm 
K 



0.5 



K300-K) = 



1.10281*10 



,"3 



2.0603*10 



~5 



1.76981-10 



k "5 



1.51851-10 
1.06784-10" 



~5 



a kg*m sec 



Mmix(x,T) = ^ 
n=0 



y(ych40,x) n n(T) n 



4 

J] y(ych40,x) m - 
m = 0 


1 + 


fKT)\°- 5 /MWJ 
lKT)J I MWJ 


0.25 


2 




/ MW \ 
8-1+ 0 

[ \ ^JJ 


0.5 



y(ych40,0.5) = 



0.00592 

0.19074 

0.76583 

0.00592 

0.03159 

0 

0 

0 



|imix(0.5 t 300 K)=0.00002*kg-m l 'sec 



ych40 =0.01183 



... F ^- „■■• : ;g f ^ ^ ... 



*0 

!. i 

OS 

w 



ru 
m 



5M0 



4.5*10 5 



M™(-5.T(fap idsp ) 4 *io" 5 
Mm K , ' Tdap idap) 



3.5M0 5 



3*10 




1100 1200 



^(• 5 ' Tda Pid*) 
^(■ 5 > Tda Pidap) 

0.998 - 




The viscosity of the CH4 - air mixtures vary 
little with composition (< 0.5%), again 
because the fuel is so dilute. We can 
therefore negfect the effect of composition by 
chosing x=0.5 in ail cases. 



7.510 2 -K 



8 lO^K 



8.5 io 2 -*: 



9 lO^K 



9.5 lO^K 



1-10*-K 



1.05- lO^K 



l.MO^K 



1.15 lO^K 



1.2 lO' K 





3.38338- 10" 5 *kgm" 


l sec~ 1 


3.53056 10" 5 -kg m" 


'■see" 1 


3.67429- 10" 5 kgm" 


'•sec" 1 


3.81491 10 ^kgm 1 


^sec" 1 


3.95267 10" 5 kg in 1 


sec" 1 


4.08779 10" '-kgm 1 


■sec" 1 


4.22047- 10" 3 kgnf 1 


-sec" 1 


4.35088 10" 5 kg in 1 


-sec" 1 


4.47917- 10" 5 -kgm~ 1 


•sec" 1 


4.60547 10" ^kg-m" 1 


-sec" 1 



WnixT(T) =exp|- 0.70216045 + 0.6408741n|^| + ' 14,2 ^ 382 K j 10 6 kg m'^sec" 1 



fimix(.5,1000 K)= 4.08779- 10 5 kg'm ^sec 1 



MmixT(1000 K)= 4.08782- 10 5 kgm *-sec 1 



Define thermal conductivity functions for light molecules 



igas r 0,l.. 4 



a3L: = 



-1.869 




0.08727 " 


- .3273 




0.09966 


0.3919 


W = 


0.09816 


-7.215 




0.08015 


7.341 _ 




-0.01013 _ 



C X: = 



1.179- 10" 4 
-(3.743- 10 s ) 
-(5.067- 10 s ) 
5.477- 10* 
1.801 10* 



dX = 





CH4 


(3.614- 10 s ) 
9.732- 10" 9 


02 
N2 


1.504 10 s 
(1.053- 10 s ) 


C02 


. -(9.1-10"*) 


H20 



1 K ' \k/ ' \k/ 



0.00001- 



W 
cm K 



I I 3 



ru 

FU 



Xmix(x,T) = ^ 
i = 0 



y(ych40,x),X(i,T) 



]T y(ych40,x).- 
j = 0 



1 + 



8- 1-h- 



MW. 
MW. 



0.5 



X(igas, 1000K) 


0.16716 kg m sec 




0.07 163 kg* m sec' 




0.06292 *kg*m sec" 


3 KT ! 


0.06788kgmsec" 


3 KT l 


0.0863 lkgmsec 


3 .KT l 



W/cm/K 



Xmix(.0, 1000 K) =0.00066 



W 

cm*K 



0.075 



0.065 - 



^( 0$ ' Td *idsp) 




0.055 



1200 



Tds l>id*> 


Xm^O^Tdsp.^) 


7.5- lO^K 




0.05362kg-msec" 3 -K: 1 


8 lO^K 




0.0562 1 kg- msec" 3 K~ ! 


8.510*K 




0.05872 kg m sec ^KT 1 


910 2 K 




0.061 nkgmsec^KT 1 


9.5- lO^K 




0.06355 kg m sec 3 -^ 1 


1 lO^K 




0.06589 kg- m-safMc 1 


1.05-10?-K 




0.06818 kg m sec ^KT 1 


l.M0*-K 




0.07043 kg m sec 3 KT 1 


1.15-10 J -K 




0.07266-kgmsec' 3 K~ l 


1.210 3 K 




0.07488 kg m sec 3 K -1 



1.01 



w ( 0 ' Td yjdJ 

^(.5,T*p idap ) 



1.005 - 



0.995 " 



0.99 




1100 1200 



The thermal conductivity varies > 1% with 
conversion. Selection of 50% conversion 
gives and error range of +.5-.5% at 1000K. 
This is likely to be smaller than the accuracy 
of measured thermal conductivity. Therefore 
we again define a composition independent 
formula to increase speed in calculating 
derived transport properties 



■es. 



WT(T) :=exp|- 7.073383 + 0.6397005- ln|^j + 65 ^ 87 K |>Wm' 1 KT 1 



Xmix(.5,1000 K) = 0.06589"m l K 1 *W 



******************** 



XmixT(1000 K) =0.06589"m l 'K 1 *W 



*»»»»*******»♦***»**»»»**»******»*♦»»» 



V|5 Determine the heat capacity and enthalpy functions of the fuel-air mixtures for the evaluation of temperature rise 

m 

hi Define heat capacities and heat of combustion: 



Cpcoef: 



19.25 .05213 

28.11 -0.000368 

31.15 -0.01357 

19.8 0.07344 

32.24 0.001924 
-9.52 0.5066 
-76.19 0.7867 
-37.18 0.8122 



0.00001197 - 0.00000001132 

0.00001746 - 0.00000001065 

0.0000268 - 0.00000001168 

-0.00005602 0.00000001715 

0.00001055 -0.000000003596 

- 0.0002729 0.00000005723 
-0.0004204 0.00000007516 J/mol/K 

- 0.0004509 0.00000009705 



Cp(igas,T):= £ Cpcoef ^ (l)" • 1 



ii = o 



Cp(0,300K)=35.66066*K 1 *J 
Cp(6,300K) = 124.01332'K" 1 *J 



Cp(5,300K) = 119.44421-K 1 *J 
Cp(7,300*K) = l68.51935-K~ 1 *J 



y(ych40,0.5) = 



c P mix(x,T) = c vW> T >y(y cM0 >% 

jk = 0 



Cpmix(0.5,998 K) = 33.66993 J K 1 



0.00592 

0.19074 

0.76583 

0.00592 

0.03159 

0 

0 

0 



Define heats of combustion 



ACpcomb(T) := (Cp(3 t T) + 2*Cp(4,T)) - (Cp(0,T) + 2 Cp( 1 ,T)) 
ACpcomb5(T) = (5Cp(3,T) + 6 Cp(4,T))- (Cp(5,T) + 8 Cp(l,T)) 
ACpcomb7(T) =(7 Cp(3,T) + 7*Cp(4,T))- (Cp(6,T)-»- 10.5Cp(l,T)) 
ACpcomb9(T) =(9 Cp(3,T) + 9 Cp(4,T)) - (Cp(7,T) + 13.5Cp(l,T)) 



AHf 



-74900- J 

0-J 

0J 
-393800 J 
-242000 J 

- 154600- J 
-119400 J 

- 103600- J 



ACpcomb(298K) = 10.38812' K 1 *J 



ACpcomb5(298-K) ^^.^^.^i 
5 6 



Define enthalpy function for individual components and AH for fuels: 

rT 



J^Hf(nn,T) =AHf + 



Cp(nn,Tx)dTx 



298-K 



CH4 + 202 ==> C02 + 2H20 
C5H12 + 802 SC02 + 6H2Q 
C7H14 + 10.502 ==> 7C02 + 7H20 
C9H18= 13.502 ==> 9C02 + 9H20 



ACpca m b7(298 K) mJ 



ACpcomb9(298.K) =8J177 . K - 1 mJ 



ACpcomb(Tx) dTx 



298K 



ACpcomb5(Tx) dTx 



^AHcomb(T):=[(AHf 3 +2-AHfJ- (AHf Q -h2-AHf 1 )]- 
AHcomb5(T) :=[[(S-AHf) 3 + 6 -AHfJ - (AHf Q + 8 -AHf J] + 

H AHcomb7(T) := [[<7-AHf) 3 +- 7-AHfJ - (AHf Q + 10.5-AHfJ] + ACpcomb7(Tx) <fTx 

J298-K 

AHoomb9(T) :=[[(9-4H0 3 +9-«HfJ- (AHf fl + 13.5-AHfJ]* 



0 



298-K 

rT 



ACpcomb9(Tx)dTx 



Hf(igas,1000K) 



-36559.44819 kg-m 2 sec 2 



22590.0387 1-kg m 2 sec 2 



21484.78794 kg-m 2 sec 2 



3.60367 K^ kg-m 2 sec' 2 



-2.15959 10 s kgm 2 sec" 2 



•»298K 



Evaluate a few values of AH for CH4, C5H12, C7H14, and C9H18: 



AHcomb(800K) =-8.00306* 10 3 J 

AHcomb7(800K) = ^ 2472 . 10 5 
7 



AHccmb5(800-K) = ^ 69221 . 10 5 mJ 
5 

AHcomb9 ( 800K )=: -6,26002- 10 5 "J 
9 



*********************** **************************** ****** •**^MHk##*#***»»*^»*##*********«**^***M*^****«*^*****^*** ******* •♦****•*•* 



Define the mixture enthalpies and determine the fuel composition needed for 51 1 F temperature rise from 1000F. 

igas=0,1..7 



7 

Hmix(x,T) = £ y(ych40,x W Hf(kgas,T) ych40 =0.01183 

kgas = 0 

Hmix(0,810.93 K) = 1.006514- 10 4 *J Hmix( 1 , 1094.82 K) = 1.0065 18* 10 4 -J conc^tr^on ^ ^ 



Hmix5(x,T) : = ^ ySCyOcSh^x^HfCkgas,!) y0c5hl2 = 0.0028968 
kgas = 0 

^ Hmix5(0,810.93K) = 1.058499" 10 4 *J Hmix5( 1 , 1094.82 K) = 1.058492' 10 4 "J 

V J Hmix7(x,T) := V y7(yOc7hl4,x) Hf(kgas,T) y0c7hl4 = 0.0021859 
T j kgas = 0 

Q Hmix7(0,810.93K) =1.075807* 10 4 -J Hmix7(l, 1094.82 K) = 1.075812-10 4 -J 

ru 



Determines C5H12 inlet 
concentration 



Determines C7H14 inlet 
concentration 



O 



Hmix9(x,T) := ^ y9(y0c9h!8,x) kgas Hf(kgas,T) 
kgas = 0 



y0c9hl8 = 0.0016878 



Determines C9H18 inlet 
concentration 



Hmix9(0, 810.93 K) = 1.083024' 10 4 *J Hmix9(l, 1094.82 K) = 1.083028* 10 4 *J 



Define momentum and thermal difusivrties: 
„I(igas,T)RT 



a(igas,T) 



Cp(igas,T)P 



P = 3.039' W Pa 
a(0,300K)=0.07813*se<f 1 -cm 2 



amix(x,T) 



, Jmix(xJ)RT 
Cpmix(x,T) P 



pmix(x,T) 



Dch4(.5, 1000-K) =0.5855*sec 1 'cm 2 
Dc7hl4(0.5, 1000 K) = O.I698' sec" 1 'cm 2 

Define Prandl, Schmidt, and Lewis numbers 



03 

W 

Pr(x,T) 



omix(x,T) 
otmix(x,T) 



omix(.5 7 1000K)=0.5352-sec 1 'cm 2 

pmix(.5,1000-K) = 1.041 17- 10" 3 gin cm 3 

i)mix(.5 > 1000K)=0.3926-sec" 1 'cm 2 

DcShlZC.S.IOOO-^^a^ST-sec" 1 cm 2 
Dc9h 1 8(5, 1000-K) =0.1 3431 sec" 1 cm 2 



Pr(. 5, 1000-K) =0.73359 



.. f 1000-K) =0.67054 Sc5(x,T) := umix ( x ' T ) 

Dch4(x,T) Dc5hl2(x,T) 



Sc5(.5,1000K) =2.00639 



Iti Sc^T) := "i**' 7 ) Sc7(.5. 1000-K) -2J1224 Sc9(x,T) 

m Dc7hl4(x,T) Dc9hl8(x,T) 



□ 



2.5 



1.5 



1 - 



0.5 



-L 



700 



800 



900 1000 



1100 1200 



Sc9(. 5, 1000-K) =2.92325 



These relations have weak 
temperature dependence and 
(presumably) very weak 
dependence on composition 



V ~ Le(x,T) 



B^xj) Le( .5 , 1000K) = 1.09404 Le5(x,T) 



_ Dc5hl2(x,T) 
omix(x,T) 



Le5(.5,1000K) =0.36563 



Le7(x,T):= Pc7hl4(X,T) Le7(.5, 1000K) =0.31726 Le9(x,T) = Pc9hl8(x,T) Le9( .5, IOOOK) =0.25095 
<anix(x,T) omix(x,T) 



1.2 



w 

m 
w 

3 

o 

Si 



^(• 5 ^id,p)o.6 



0.4 - 



0.2 



_L 



700 800 900 1000 1100 1200 



5= l ************************** 



ru 

nJ Can remove compositional dependence for most property fundctions 

□ 
M 1 



pmixT(T) =MWmix(0.5) 



RT 



i>mixT(T) - 



HmixT(T) 
pmixT(T) 



The Lewis numbers show 
approximately the temperature 
rise at the wall of an active 
catalyst relative to the adiabatic 
temperature rise at that point. 



Clearly temperature (and composition) 
have little influence on the AT at the 
wall, while the nature of the fuel 
molecule has a strong effect (from 
0.25* ATadiab for C9 to 1 .1 *ATadiab 
forCH4) 



#**#**•*•♦**********«*##*###**#♦**#***** 



ScT(T) 



pmixT(T) 



Dch4T(T) pmix(0.5,T) 



Sc7T(T) := 



HmixT(T) 



Dc7hl4T(T)pmix(0.5,T) 



These functions speed calculations of Nu numbers 
and heat and mass transfer coefficients 



This ends the extended section that defines gas properties for the combustion of several fuel air mixtures. 



Define honeycomb configurations 

Diao = 9in Diao = 22.86 -cm Diai -4 in Diai = I0.l6-cm Lhc=3.5*m Lhc=8.89'cm 

With bypass change Oiao to 8.5 in Bypass = 0 Bypass - 0, 0.2. 0.4 Lhc - 3.5, 4 in 



Honeycombs: Represent as square channel metal alloy wall 0.05 mm, uncoated 
Cells 1 := 100 in" 2 dwalil =0.009 cm Cells2 = 200 m 2 dwaJ12 = 0.009- cm 

dcelll =(CeIlsl)" 05 -dwaUl dcell2 = (Cells2)~ 0 5 - dwaI12 

dcelll =0.245 -cm ksl =0.002- cm dcell2 = 0.17061 -cm ks2 =0.002- cm 

ocelli = dceM 2 Cellsl eceDl =0.93039 ecell2 = dcell2 2 -Cells2 eoeI12 = 0.90229 



Describe 100 and 200 cpi 
dimensions, although not 
expected to be used 



Cells4 =400 in 2 dwa!14 =0.007 cm 



dcell4 :=(Cells4)" 0 5 - dwaH4 dcell4 =0.1 2 'cm 
dcell4 =0.12 'cm ks4 =0.001 cm 



p. ece!14 =dcel!4 2 Cells4 £ce!14 = 0.8928 



~J3 



Cells6 =600 in" 



dwa!16 =0.007 cm 



dcell6 =(Cells6)'° 3 - dwall6 dcell6 = 0.0967 *cm 
dcell6 =0.0967 -cm ks6 =0.001 cm 

£&eU6 = dcell6 2 Ceils6 eee!16 = 0.86955 



Packed Beds: Cylindrical (conical) geometry 



UJ Dim '= 6- in 

L, Dim* = 15^4 'cm 

Q 
DO 

LrT: Smono ■ = --(Diao 2 - Diai 2 ) 

* y 4 

p 



Doutr =7 in, Wbed =7 in 

Doutr = 17.78 'cm Wbed = 17.78 -cm 



Sbed:= ,.|^5Lt5^j.Wbed 



Lbed =0.75- in 



Smono 
Sbed 



=0.35714 



Define initial combustion conditions 



Adjust U0 to fit U(Tin)-7 m/sec w/o bypass flow 



Uref =7- 



m 

sec 



Tin =(537.778 + 273.15)K Tin =810.928'K 



M0 = Uref- pniix(0, Tin)- 



UO 



. MO (1 - Bypass) 



f((9m) 2 -(4in) 2 ) 
4 



1 



P= 3.039* KT -kg-m '*sec 



MO =0.29601*kg"sec 



Tout =(821.667 + 273.15)K 
Tout = 1094.817-K 



U(T)=U0--^ 



Ubed(T) 



M0 



pmix(0,Tin)Sbed 



U(Tin) = 7*nVsec 



Ubed(Tin) =2.5*m*sec 



G0=U(Tm)pmix(0,Tm) 



GO = 8.98749*kg-ni -sec 



V 



Evaluate Reynolds numbers for use in calculating the momentum, heat, and mass transfer functions; 



Rehc4{T) 



OOdoelU 
eccI14mmxT(T) 



Rehc4(Tin) =339.11451 



Rehc6(T) 



G0-dceU6 
ecel!6mnixT(T) 



Rehc6(Tm) =280.56407 



For laminar flow with entrance effect (adapted from Kays & Crawford, "Convective Heat and Mass Transfer"): 



cfiric4(T,z) 

Reho4(T) 



1 + 0.046263- Rehc4(T) 



dce!14 



10.45363 



14.227 



Refae6(T) 



1 + 0.046263 -|Rehc6(T)~^J 



10.45363 



^ cfric6(T,z):= 

yp 

m 

! . 3 
! : * 

Calculate monolith pressure drop in several monoliths with entrance effect: 

m 



Average coefficient of friction 



2 idap r ^.5 cnn- 1 idsp cm 



Udsp ;dap : = (4 + idsp) m sec" 



00 

ru 

w 



0.08 »" 

c^lOOfrK.z^ 

14.227 
Rd>c6( 1000-K) 

14.227 o.06 - 

Rchc4( 1000-K) 



0.04 




0.06 0.08 



z idap 



The entrance effect for AP of a channel in laminar 
flow is significant even to 3.5 in. length. 



v Tavg- 



Tin+ Tout 



cficx(x,z) :=li??Z .(i + 0.0462363 (xz))° 45363 



D £ _dceU6 Uds Pid» 
Rex6. = = — 

cceU6 umixT(Tavg) 



Rex4 



tdsp 



. dcelM Uds Pj,hp 
ecell4 umixT(Tavg) 



\ Lhc / 2 \ ecelM 

\ ^ Lho / 2 \ eceD6 / 



dcelM 



\ ^ Lho I 2 \ eceI16 / dcell6 



0.006 



^li- 
fe* ! 



ri 




*************************** 



APhc6(T) = 4cfric6(T,Lhc) 



pmixT(T)(^>] 2 

\<fceH6/ Lhc 



AP/P < 0.5% to about 1 1 m/sec 



#♦**##♦♦♦♦♦*♦#♦**•****•**♦«♦♦♦♦*♦**#*#*♦*##***♦♦#*♦♦##****♦*#* *****t 



dce!16 



T«vg=952.8725*K 



U(Tavg) =8.22528"msec 



^(Im) =0 00271 



APho6(Tout) 



= 0.00441 



APfac6(Tavg) =0()0352 
P 



APhc4(T) =4cfric4(T,Lhc) 



pmixT(T)-(W) 2 

\eceU4/ Die 



dceU4 



APho4(T«vg) =0Q0228 



******* ******** ********************* ************************ ******* 1 



>»»»*»*****»»*♦*♦** *1 



r************ **** i 



Find 1-0 heat and mass transfer rates via Nussett formulas 



Averaged Nu for 600 cpi square channels (fully developed 
laminar flow) taken from Rosner: 



NuhavgR6(z) =2.976 



1 + 



7.6z 



\Reho6( Tavg) Pr( .5 , Tavg) dcell6 



7.6(2cm) 



Refae6(Tavg) Pr(.5,Tavg) dcel!6 



NuhavgR6(Lhc) =2.98649 



= 0.85369 



NumavgR6(z) =2.976 



l 

818 



1+1 



7.6z 



\Reho6(Tavg) ScT(Tavg) dce!16 



NnmavgR6(Lhc) =2.98441 



Km 



; Num5avgR6(z) = 2.976- 



1+1 



7.6-z 



\Rehc6(Tavg)- Sc5( .5 , Tavg)dcell6 



)1 



w 



a 



* Num7avgR6(z) =2.976 



1 + 



( — ) 

\Rehc6(Tavg)Sc7T(Tavg)dceU6/ 



Num5avgR6(Lhc) =3.11 102 



Num7avgR6(Lhc) =3.16167 



b 



Nam7avgR4(z) = 2.976- 



1 + / 2^5 ) 3 Num7avgR4(Lhc) = 3.3902 

\Rdic4(Tavg)Sc7T(Tavg)dcell4/ J 



Num9avgR6(z) =2.976 



1 + 



7.6-z 



\Reho6(Tavg) Sc9( .5 , Tavg) dcel!6 



Num9avgR6(Lhc) =3.2775 



The Nusselt numbers for mass transfer show small (< 5%) compensation for the low diffusivity of the higher hydrocarbons. 



Use exponential functions implied by the definition of average (integrated) Num functions to calculate exit mole fractions. 



yc7h!4<T,z) ::y0c7hl4.exp( ^Nmp7avgR6(z) 

\Rehc6(T) So7T(T) dceU6 



yc7h!4(Tavg,Lhc) =() 13?()5 
y0c7h!4 



ych4(T,z) :=yeh40exp( ' ^ ^Nm»uvgR6( 2 ) 

\Rehc6(T)ScT(T)-dcelI6 



ych4(Tavg,Lhc) mQMim 
ych40 



yc5hl2(T,z) ^y0c5hl2expf 4 ^^av^z) \ 
\Relio6(T).Sc5(0.5,T).dcell6/ 



yc5h!2(Tavg,Lhc) 
y0c5h!2 



= 0.10114 



09 



yc9h!8(T,z) =y0c9h!8 exp 



o.oi 



yc5hl2(T«vg,z id ^ 



0.001 



1*10 



-4zNmn9avgR6(z) \ 
Rfihc6(T)Sc9(0.5,T>dcell6i 



J. 



1 



0.02 



0.04 0.06 



yc9h!8(Tavg,IJic) =()1%04 
y0c9h!8 




0.08 



0.1 



z 9 =0.095-iii 



z 9 =3.74016 in 



Summary if 600 cpi Metal Monolith Performance 

%Bypass Length(in) Diao(in) Yout/Yin C7H14 Yout/Yin CH4 



0 
0 
20 
20 
40 
40 



3.5 
4 

3.5 
4 

3.5 
4 



9.0 
9.0 
8.5 
8.5 
8.5 
8.5 



0.13705 
0.10679 
0.11888 
0.09055 
0.06266 
0.04317 



0.001533 
0.000611 
0.000906 
0.000335 
0.000089 
0.000023 



AP/P 

0.00352 
0.00400 
0.00324 
0.00368 
0.00240 
0.00273 



* Bypass = 0 

Lhc =3.5 'in 
Diao =9 in 

These values are conservative 
and do not account for the 
higher average temperature 
with the higher F/A ratios and 
do not account for 
hydrothermal entry. 
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0.0 1 0.02 0.03 0.04 



zx =25.25*cm 



APhc6(Tavg)>zx =0010()0 
PLhc 



y ch4(Tavg >Z x) =106492<10 - 
ych40 



PLhc 



y0c7hl4 



y^hlSCTavg^zx) =14421 g. 10 - 
yOc9hl8 



Constant Flow 8.2 m/sec 
Variable Length to 89 cm 

U(Tavg) = 8.22528'nrsec~ 1 
lOLhc =88.9'cm 



Conversions at AP/P=0.01 
CH4 = 1.06*10 A -8 
C7H14=4.8*10 A -3 
C9H1 8=1 .44*1 0 A -2 
z=25.25 cm 



Conversions at AP/P=0.02 
CH4 = 1.14*10M6 
C7H14=2.41*10 A -5 
C9H1 8=2.21*1 0M 
z=50.5 cm 



The 600 cpi monolith has low pressure drop at 8.89 cm length. If it's length were extended to give 1 .0% and 2.0% 
fractional pressure drop, then C7H14 conversions would increase to 99.5% and 99.998% respectively 



\ ^ Set up properties for CO to calculate the conversion of dilute CO and hydrocarbons downstream of the turbine 



DcoI(x,T) 



P 
Pref 



Dco7(x,T) 



E 



y(ych40,x). 



jgas = 0 



Diff(8jgas,T) 



P 
Pref 



£ y7( y 0c7hl4 > x) jgM 



jgas=0 



Diff(8,jgas,T) 



Dcol(.5,Tin) = 1.12854'sec 1 'cm 2 Pref =l'atm 

Dch4(.5,Tui)=0.41l95sec" 1 "cm 2 P=3atm 

Dco7(.5,Tin) = 1.12396'sec _1 -cm 2 Pref=iatm 
Dc7hl4(.5,Tin) =0.1 1846 -sec 1 'cm 2 P=3*atm 



0.0003 



W - 



*01 

en 



Deolf.J.TApy^ 

Dco7(.S.Tdqi ida|1 ).0001 




ru 
m 

o 



DcoT(T) =exp 



10.800864- 1.6357427 



157906 



K 



10' 4 m 2 ! 







Dcol(.5,Tdsp idap ) 


750-K 




9.90331 10" 5 m 2 sec" 1 


800-K 


1.10324 10" 4 m 2 sec" 1 


850K 
900 K 
950K 


1.2208 10" 4 m 2 sec 1 


1.34291 10" 4 m 2 sec 1 


1000K 


1.46949 10 4 m 2 sec 1 


1050K 


1.60044 10* 4 m 2 sec" 1 


HOOK 


1.73572- 10" 4 -m 2 -sec" 1 


1I50K 


1.87523 10 4 m^sec" 1 


1200K 


2.01893 I0" 4 m 2 sec 1 






2.16675- 10" 4 m 2 sec 1 



Dcol(.5,Tavg) = 1.47689-sec 1 *cm 2 
DcoT(Tavg) = 1.47689'sec' 1 -cm 2 



dreco "•= 17 - in 
dreci ~6 m 

M0 = 0.29601 -kg" 
M0 



GrecO 



Arecup 



Prec 1 atm 



dreco -43.18 'cm 
dreci = 15.24 'cm 



Arecup 



dreco 2 - dreci 2 



GrecO = 2.30904*kg*m~ 2 'sec 1 
MWmix(0.5)Prec 



pmixrec<T) :=- 



RT 



Lrec -3.5 in 
Lrec =8.89*cm 

Arecup =0.1 282* m 2 



GO =8.98749 'kg' m" 2 ' sec 1 



pmixrec(Tavg) =0.36422'kg'm 



Urec(T) :=- 



GrecO 
pmixrec(T) 



Rerec4(T) := 



Rerec6(T) 



GrecO- doei!4 
ymixT(T) 

GrecO dce!16 
jimixT(T) 



Urec(Tavg) = 6.3 3964 'or sec 1 

Rerec4(Tavg) =69.96214 
Rerec6(Tavg) =56.37495 



SccoT(T) :=- 



WnixT(T) 



DcoT(T)-pmixrec(T) 

Schmidt numbers are 
independent of pressure 



Tavg=952.8725'K 



cfiicrec4<T,z) 



. 14.227 
Rerec4(T) 



1 + 0.046263- Rjerec4(T) 



dceD4\ 



0.45363 



cfncrec6(T,z) - 



14.227 
Rerec6(T) 



<kell6 



1 + 0.046263 |Rtoc6(T) 



10.45363 



*8 

^ Calculate monolith pressure drop in several monoliths with entrance effect: 



i . a 

w 



pmixrec(T) 



APrec4(T) = 4cfricrec4<T,Lrec) 



/ Urec(T) \ 2 

\ ecelW / Lrec 



pmixrec(T) 



03 APrec6(T) =4 cfricrec6(T,Lrec) 



2 dceU4 

/Urec(T)\ 2 
\ ccell6 / Lrec 



dceU6 



APreo4(Tay g> =0.00557 
Prec 



*B*Sl2!ti. =0.00898 

Prec 



Numavgrec6(z) = 2.976- 



1 + 



7.6z 



Rerec6( Tavg)*ScT( Tavg)*dcell6 



Numavgrec6(Lrec) =2.97616 



Numavgrec4(z) =2.976 



1 + 



7.6z 



Rerec4(Tavg) ScT(Tavg) dcell4 



Numavgrec4(Lrec) =2.97649 



Numavgrec7R6(z) =2.976 



l 

8l» 



1 + 



7.6-z 



Rerec6( Tavg ) Sc7 T( Tav] 



gHcelte) 



Numavgrec7R6(Lrec) =2.98024 



Nnmavgrec7R4(z) =2.976 



1 + 1 



7.6z 



\Rerec4(Tavg)Sc7T(Tavg)dcell4 



Numavgrec7R6(Lrec) =2.98024 



NumavgreccoR6(z) =2.976 



1 + 



7.6z 



\Rerec6(Tavg) SocoT(Tavg) dccll6/ 



NumavgR6(Lrec) =2.98441 



NumavgreccoR4(z) =2.976 



l 

8*1* 



1 + 



7.6z 



\Rerec4(Tavg)SccoT(Tavg)dcell4/ 



NumavgreccoR4( Lrec ) =2.97663 



The average (integral) Nusselt numbers at 1-atm with the larger flwo cross section show very little contribution of the therma 
entry length, i.e. the flow velocities and heat and mass transfer coefficients have values close to those in long channels. 



p Now determine CH4, C7H14, and CO conversions in 400 and 600 cpi honeycombs, assuming 80% C7H14 conversion 
.fj downstream of the combustion monolith and upstream of the turbine 



2 ych4rec6(T,2) :=ycMCT«v > ,Uio)-«p( ' 4 •*M""'g«*(*> | 
y K \Rerec6(T)ScT(T)dcell6/ 

^ ycfa4nxtf(Tin,Lrec) =0 ych4(Tavg,Lhc) = 1.81356' 10 5 



y ch4rec6(Tavg,Lrec) = 2-733018 . l0 ~n 
ych4(Tavg,Lhc) 



Bypass = 0 



Lhc=3.5 'in 



Diao =9 'in 



m yc7hl4n»6<T,z) ^O^ycTh^Tav^LhQ cxpf 4 ) ycThUn^Tin,!^) . 

■ ^ ^ * J y ^ * ^\Rerec6(T) Sc7T(T) <toell6/ 0.2 yc7h!4(Tavg,Lhc) 



f=t 



yc7hl4rec6(Tm,Ln») =3.42819*10 



0.2-yc7hl4(Tavg,Lhc) =5.99172*10 5 



ych4rec4(Tin,Lrec) =8.84152*10 



13 



yc7hl4rec4(T,z) =0.2 yc7h!4(Tavg,Lhc) exp 



/ -4 Z'Nmnavgrec7R4(z) 



yc7h!4rec4( Tin, Lrec) =Q (K)773 
0.2yc7hl4(Tavg,Lhe) 



\Rcreo4(T)-Sc7T(T)dcell4^ 

^ yc7hl4rec4(Tin,Lrec) =4.63263* lo" 7 0.2yc7hl4(Tavg,Lhc) =5.99172* 10~ 5 0.4 7 yc7hl4(Tavg,Lhc) -8.3884*10~ 4 



yeorec^z) ==M-7-yo7hl4(T W , > LhD)-« > ( 4 2 ^g^K^z) [ 
* V ^ e ^\Rerec4(T)SccoT(T).dceU4/ 

J ^\Rerec6(T)SccoT(T)dceU6/ 



ycareo4(Tavg,Lrec) = 3.06761- 10 



n 



ycorec6(Tavg,Lrec) = 2.9626* 10 



15 



zx. idap -(idsp + 0.05) cm 



0.001 



Q 



u 

m 



HI 

Q 



ycorec^Tavg.z*^ 
yc7hl4roc4(T«vg,2X idsjj ) J 




Sc7T(Tavg) =2.31842 

SccoT(Tavg) =0.73627 
ScT(Tavg) =0.67113 



ycorec4( Tavg, Lrec ) 
(0.4-7 yc7h!4(Tavg,Lhc)) 



= 3.65696-10 



yoorcc6(Tavg t Lrec) ^ 5317g . 
(0.4-7yc7hl4<Tavg,IJic)) 



10 



12 



This analysis assumes 80% homogeneous conversion of fuel exiting 
a 3.5in. 600 cpi monolith without bypass flow and assuming 40% 
production of CO in the post combustion monolith region. 

CO conversion in a 400 cpi recuperator monolith is >99.99%, 
while cyclohexane conversion is 99.2%. 

The final CO and C7H14 levels are < 1 ppb and .5 ppm while post 
turbine concentrations are taken as 840 and 60 ppm respectively. 
CH4 conversion also should be very high and use of a 600 cpi or 
longer monolith would give higher HC conversions (x20 for 600 cpi). 

Post turbine pre recuperation monolith pressure drops are modest 0.56% 
for 400 cpi and 0.90% for 600 cpi monoliths. 



********** ************ **«**i 



*********** ************** 



******* *****«*«#»«***#*♦*♦♦♦***»**#*#•**♦*♦*****♦♦*♦#**##*### **#^** * * ************ *** • *************** 



v - An alternative method of calculating performance uses local Nusselt numbers for heat and mass transfer to determine the 
wall temperature and combustion rates at the wall. First we define the local Nu function then check conversion assuming 
feist kinetics and compare with the previous averaged (integrated) Nu formula. 

Local Nu(x) taken from Kays & Crawford for hydro & thermal entrance laminar with circular channel: 
Numono(z) =(1 +exp(- 1.09689 - 0.38 14021n(z) - 8.287774-z 0 5 )) 3. 6<Nnmoao(0. 1 152163) =3.82723 

^ _ Rehc6(Tavg) -Pr(.S,Tavg)-doeU6 

ZZCMn> " 2 zzeh*6 = 8.90257 -cm =0.99859 

zzchai6 T 

Rehc6(Tavg) =252.34858 
itest : = 0 50 dce!16 = 0.0967 4 cm 

ztest_ + dwal!6 

ztesT ^ = ftest*0.05cm + 0.05*cm zz._ : : 



itest 



zzchar6 




. „ * tJ . Rchc6(Tavg) ScT(Tavg) dcell6 
Rehc6(Tavg) = 252.34858 doeU6 = 0.0967 'cm zz6ch4 = ~ &/ ■ — 

Pr(.5,Tavg) =0.72969 

Rehc6(Tavg)Pr(.5,Tavg)dcell6 = 17.80513 -cm 

z26c7hM . Reho6(Tav g ) Sc7T(T>v g ) dceU6 



XmixT( T) Numono 



2 

z 



. /T , , \zzchar6 

b(T,z) 

doell6 

Dch4T( T) Numoaof — - — J 

I , kgpb(T,z) := \zz6ch4/ kgpb( Tin, Lhc) = 15.5934 sec 1 "cm 

dcel!6 



• — — — m 



Dc7hl4T(T)Numono 



kgpb7(T,z) =- 



\zz6c7hl4/ 



dceU6 



kgpb7(Tin f Lhc)=4.50626 sec 1 'cm 



Prater Temperature for CH4 at the 
inlet with fully developed flow 



kgpKT«n > Lhc)^ ych40 -^ C ^ b(Tin) =309.73474'K 
RTin h( Tin, Lhc) 



Prater Temperature for C7H14 at 
the inlet with fully developed flow 



kgpb7(Tm,UK ; )-g- y 0c7hl4 - AHcomb7(Tk > =90.38489'K 
eK RTin h(Tin,10cm) 



***************************************************************************** 



Determine maximum conversion of fuels with very fast catalytic combustion rates (zero fuel concentration at the wall) 

Tin+ Tout 



iz =0.. 100 



ych46 ft : - ych40 



yc7hl46 = y0c7hl4 



Tavg : = - 



W 

m 
W 

S 

o 

I i 5 

ru 

m 



^ . llr . Lhc 7 

z6. - dwalK) + iz 

u 100 7 



ych46. 



12+1 



1 - kgpb(Tavg,z6. ) ych46 ; , 

*-M *• a/ l00 7 dceU6 U(Tavg)/ u 



Lhc— =8.89 'cm 
7 



3.5 & 4-in. monolith length 



yCh46 - =1.2447 10- 



ych40 



ych4^Tavg,Lhc*- 
ych40 



= 0.00153 



^146^,:= II 



kgpb7(Tav g ,z6. )i* !-* ^-) yc7hl46. 

™ V 6 ^ 100 7 dce!16 U(Tavg)/ 



yc7hl46 



y0c7hl4 



i2i =0.06145 



•*♦*«♦**♦♦****#**♦**♦♦♦♦**************♦**♦♦* 

By-pass length yC7H14ow7yin yCH4outfyin 



0 


3.5 in 


0.06145 


1.2447*10 * 


0 


4 in 


0.04256 


3.288 10" 5 


0.2 


3.5*in 


0.0502 


5.96324* 10" 5 


0.2 


4 in 


0.03374 


1.3974 10* 3 


0.4 


3.5*in 


0.02054 


2. 17587- 10' 6 


0.4 


3.5 in 


0.01207 


2.92928- 10" 7 



Tavg=952.8725*K 
Lhc =3.5 -in 
Diao = 9 *in 
Bypass = 0 

Nu function with circular channel hydothermal entry length 
shows 93.9% to 98.8% conversion of C7H14 



****** ****** *****#*************** ************* 



Pressure Drop Calcuations for a Packed Bed 
given: e - void fraction; 

dp - effective particle diameter 
U - mean superficial veloci ty 

vHcinematic viscocrty fbed(x) :=— + 1.75 

p- density x x - *>&<^»*~ 



I IMF 

At T - 900K and 10 atm pressure, 



pmix(0.0,300K) =3. 47058- kg- m 
pmix(.5,Tavg) = 1. 09267* kg- vP 

j = 0..3 



umix(.5,Tavg) = 0.00004 'm 2 -sec 1 



XT. 



hi 



. _ 125 . 
dp. in 

> j+1 
e=0.5 



i =1.. 10 
Utx. = im-sec~ 



Lbed=0.75'in 
oO -omix(.5,Tavg) 



pO = pmix(.5,Tavg) 



0.3175 
0.15875 
0.10583 
0.07938 



'cm 



Rebed :=Utx - 



>o 



1 1)0 (1 - E) 



□ 
S3 

ru 



0 



AP. . =Lbed fbed (Rebed V 

io v »o/ 



1-e 



p0-(Utx.) 2 



Utx = 



j3 
[4 
|5 
[6 
[7 
[8 
[ 9 
10 



-1 



So 


0 


0 


o 1 [Ho 


0 


0 


0 


|| 175.19103 


87.59551 


58.39701 


43.79776 \ IS 0.00067 


0.00179 


0.00335 


0.00536 j 


|g350.38206 
ff525.57309 


175.19103 
262.78654 


K lJ^7?4p2 
175.19103 


.?L5?55lJ H 000226 
131.39327! M 0.00474 


0 I 0054_ 
0.01081 


0.01822 


0 :9]*Hj 
0.02695 ! 


|| 700.76412 


350.38206 


233.58804 


175.19103! IS 0.00814 


0.01804 


0.02972 


0.043181 


K 875.95515 


437.97757 


291.98505 


218.9887<Ap=|H 0.01243 


0.02708 


0.04395 


0.06304! 


|| 1051.14617 


525.57309 


350.38206 


262.78654 


M 0.01764 


0.03794 


0.0609 


0.08651! 


|| 1226.3372 


613.1686 


408.77907 


306.5843 


|| 0.02375 


0.0506 


0.08056 


0.113621 


H 1401.52823 


700.76412 


467.17608 


350.38206 


|| 0.03077 


0.06508 


0.10294 


0.14435! 


11 1576.71926 


788.35963 


525.57309 


394.17982 


11 0.03869 


0.08137 


0.12804 


0.1787 j 


S.1751.91029_ 


_875.95515_ 


_583.9701_ 


_437.97757_ 


H 004752 


0.09947, 


0.15586_ 


0.2 1668 j 



'atm 



0.06 



0.05 



P 



0.04 



— P 0.02 



0.01 
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UtX; 



SAP(e,T,z)=zfbed 



Ubed(T)dp 



[ OO(I-E) 



1 - s pO (Ubed(T)r 
2 dP ; 



£=0.45 



dp = 



ru 
ru 



0.3175 
0.15875 
0.10583 
0.07938 



AP(.45,Tavg,Lbed) AP(J,Tavg,Lbed) 
P P 



"cm 



0.00173 



0.00407 



0.00703 



0.0106 



0.00113 



0.00263 



0.00449 



0.00673 



10 



dp =0.04167 'in [/ , 
dp, =0.0625 'in J^' 



dp =0.125 "in 



Lbed=0.75'in 



Ubed(Tavg) = 2.5'nrsec 



For larger particles use axial beds with greater length 

k =0 dp n =0.125 *m 



APaxial(e,T,z) -z fbed 



U(T)dp k 



oO(l-c) 



1 - £ pO(U(T)) z 



v 

^ — 



U(Tavg)=8.22528'm"sec 1 
APaxiaKfi.Tavg.l-in) =0.0657 'atxn 



P=3*atm 



APaxial(£,Tavg,l m) = QQ2 J 9 1 /8th jnch gjye gbout 2% pressure drop jnch 



tt*ttmM*mtHt«n**tmt*ttttt«*tt***m«t***tmmt«HM**ttm«t*m******mtt*«tmt*»*ttt*ttttttM*t** 



Determine Nu values for heat and mass transfer between particles and the gas phase in the bed 



Rcpartd) :s -J*5-GQ4L_ 

umixT(T) (l - e) 



Repart(Tavg) = 



398.16156 
199.08078 
132.72052 
99.54039 



dp = 



0.3175 
0.15875 
0.10583 
0.07938 



cm £=0.45 



Pr(.5,Tavg) =0.72969 



k=0 



Sc(.5,Tavg) =0.67112 



NnhpartR(T) =1.34 



Pr(.5,Ty 



k 0.4 



0.4(Repart(T) k ) 05 + 0.2(Repart(T) k ) 3 



NuhpartR(Tavg) =22.21552 



havgpart( T) 



o 
w 

i UJ Effective boundary 
lip layer thickness: 

w 

E 



. XmixT(T)Nnlq)artR(T) 



havgpart(Tavg) =445.59219 m" 2 K 1 'W 



NuhpartR(Tavg) 
dp k =0.3175 -cm 



= 142.9181pm 



ft********************* 



m 

m NmnpartR(T) := l.34-ScT(T) 04 

1 «? 



0.4 (RepaittT^ 2 + 0.2(Repart(T) k ) 



dce!16 



NuhavgR6(Lhc) 



doeU6= 0.0967 -cm 



NumpartR(Tavg) =21.48438 



= 323.77508 'Mm 



kgpartd) .^(TVNmnpartRd) 
*<*P k 



kgpart(Tm)=72.76151-fiec 1 'cm 



Nupait7R(T) := 1.34Sc7T(T)' 



0.4 



0.4 (RepartCDj 0 5 +■ 0.2- (Repartd)^ 3 



Nupart7R(Tavg) =35.27582 



wut,^ _Dc7hl4T(T)-Niroait7R(T) i^t, t - x -*a aa**<. "> . 

kgfl>7( T) — — — — - kgfb7( Tin) = 34.44665"sec 'cm 

*<*P k 



1 - kgpart(Tavg) 



Lbed 6 (1- e) 1 



50 dp fc Ubed(Tavg) 



ych46. 



ych46 



31 



ycb46 0 



=0.02369771 



ych4p(T,z) -ych40exp 



6 (1 - c) Num P artR ^ T ) z 
eRepart(T) k *ScT(T) dp k 



ycMp(T,vg,Lbed) =Q(?2908261 
ych40 



yc7hl4(T,z) -y0c7hl4exp 



.6(l-e) Nnp»rf7R(T) _z 



£Rq>ait(T) k So7T(T) dp k 



yc7W4(TavgJ^ =018611 
yOc7h!4 




dpart(cm) yout/yinC7H14 
0.3175 0.18611 
0.15875 0.0116 



youtfyinCH4 AP/P 
0.02908261 0.00113 
0.0000846 0.00263 




dp k =0.3175 -cm 

Ubed(Tavg) =2.5'm*sec" 
£=0.45 

Wbed = 17.78 *cm 
Lbed = 1.905 -cm 



i . i 

W 



in 

fcsz? 



Determine hydrocarbon conversions in the packed beds vs pressure to estimate effective conversion 
at equivalent pressure drop 



yc7hl4dp2(T,z) =y0c7hl4exp 



-6(1 -e). 



1.34Sc7T(Ty 



0.4 



0.4(Rep«t(T) 2 ) 2 + 0.2-(Repart(T) 2 ) 3 



E*Rq>art(T) y Sc7T(T) 



z 



ych4dp2(T,z) =ych40exp 



0.4 



6(1 -e) 



I.34ScT(T) u 



0.4(Rq>art(T) 2 ) 2 + 02(R^)ait(T) 2 ) 3 



e-Repart(T) ScT(T) 



dp 2 



APdp2(Tavg,z) = z fbed(Repart(Tavg) 2 ) 



1 - e pO(Ubed(Tavg)r 
e 3 " 



Ubed(Tavg) =2. 5 m* sec 



1 

0.1 
0.01 
0.001 

yc7hl4dp2(T«v 8 ,z idaD ) _^ 
y0c7hl4 

y cMa P 2 ( Tgv g>W ri0 -5 

ycb40 

1M0~ 6 



1*10 ' h 
IM0~ 



IMo"" 9 







1 




i i 
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0 




0.005 




" 0.01 0.01S 0.02 



W 

a 



m 
o 
u 



y0c7hl4 

ycMO 



P 




dp 2 = 0.10583 -cm 
dp 2 =0.04167 in 



z 9 = 3.74016 in 



Lbed=0.75*in 



dp 0 =0.3175"cm 
dp Q =0.125 *in 



z 9 = 3.74016 -in 



Lbed =0.75 -in 



Ubed(Tavg) = 2.5 'or sec 



For oblique flow through cone-shaped beds, large particles give a very low AP (even with thick beds) and high conversion 



zxx - 2.71 cm 

APdp2(Tavg,zxx) =Qm{m 
P 



yc7hl4dp2(Tavg,zxx) „ 1J2434 . 10 ~5 ycb4#2(Tavg,zxx) m $ 4$19y lQ ~n 
yOc7hl4 ych40 



For AP/P = 1% t C7H14 conversion = 99.9987% with fine particles, dp = 1.06 mm 
zy "2.0- cm 



Repaxial(T) 



U(T)dp n 



umixT(T) (1-e) 



Numpaxial(T) = 1.34 ScT(T)' 



0.4 



0.4 (Repaxial(T)r +- 0.2 (RcpaxiaKT)) 



r .i 2i 



6- ( I- E ) NpaygfayT) z_ 

E-Rep«xial(T) ScT(T) <fo 0 



Nnmp7axial(T) = l.34Sc7T(T) 0 - 4 [o.4(Rqjaxial(T)) 2 + 0^(Rq)axial(T)) 3 ] 
ych4pixi»l(T,z) =ych40-exp|- 



yc7hl4p«ul(T, Z ) =y0c7hl4 exp[-6 (l - £) Nnmp7axial(T) z_ 

' *Rep«orial(T)Sc7T(T) dp Q 



it 1 

hi 



□ 

m 
m 
ry 

! t 

WHIT 



0.1 - 



yc7hl4pmal(Tavg,2 i<|ap ) 001 
y0c7M4 

y cM P gdrf ( Tav 8' 2 id«p) 

ydi40 0.001 



1*10 



"4 _ 



1*10 3 




0.02 



0.04 0.06 
AP«rid(s,Tw 8 ,z id J 
P 



0.08 



Pressure drop increases 
approximately with 
thesquare of flow velocity 
while conversion 
decreases exponentially 



U(Tavg) =8.22528-nrsecf 1 zy =2.32 ran 

APaxj^Tavg,;?) = qq2000 yeh4paxud(Tavg,zy) =006g88284 yc7hl4paxial(T«vg,zy) _ 7RmR 

y ch40 — y^ 028038 

For flow perpendicular to the fixed bed, even with large catalyst particles, the equivalent pressure drop for 
/ high C7H14 conversion is very high. Even with AP/P = 2%, the C7H14 conversion is only 72%. For high conversion with 
practical fixed beds, the superficial bed flow rate must lower than about 3 m/sec to give high conversion with low AP 



Define Catalyst Functions: 



nsph(*) — • 



l 



♦ \tnh(*) ♦ 



TOL =0.0001 



R=82.07601-K 1 *atmcm 3 atm*cm3/mol*K 
Rem =83 143000-gmcm 2 sec ~ 2 KT 1 



nfpK*) 



= tanh(*) 
♦ 



Physical Properties for Pd/AJ203 Washcoated cCatalyst 



Scat :=85-m 2 -gnf 1 
Vnpore -0.15 cm 3 gm" ! 



pcatpart = 2.5 gm cm 3 



Loat =0.0020cm 
j« Washcoat thickness - 20 microns 



P.. ..^ Vnpore 
|jj dppore 6- 



Scat 



m 

lU ***** 



efipore =0.375 
xnpare =3.0 



dpcat =0.0002 cm 

Size of Y-AI203 particle within washcoat - 2 microns 



dppore = 1.05882* 10 6 'cm 



□ 



DeflCH4(x,T) 



ejipore 



1 



Dch4(x,T) 



DeffC7H14(x,T) := 



djipore- 
tjipore 



-RT 



1 



Dc7hl4(x,T) 



djipore- 



dppore- 



-RTin 



= 0.10953'sec 1 *cm 2 



Dch4(.5,Tm)=0.41195*sec 1 'cm 2 



DeflCH4(.5,Tin) =0.004 19' sec 1 cm 2 



dupore- 



-RTin 



= 0.04428- sec 1 *cm 2 



Dc7hl4(.5,Tin)=0.11846*sec 1 *cm 2 



DefiTC7H14(0.5, Tin) =0.00 164- sec 1 W 



DmacCH4(x,T) ™ Dch4(x,T) 
DmacC7H14(x,T) =^ Dc7hl4(x,T) 
kOT =0.01 gm" ' sec 1 Test value for kOT 



DmacCH4{.5,80aK)=0.06442*sec~ l 'cm 2 
DmacC7H14(.5 t 800K)=0.0I85rsec~ l W 
Macropore porosity - 0.4 and tortuousrty - 2.5 



Microscopic Thiele Modulus for CH4 



,x,T) := J 



pcatpart 



kOT 



DeffCH4(x,T) 



dpcat 
6 



♦Kk0T,.5, Tin) =0.021 



Macroscopic Thiele Modulus for CH4 

<hnac(kOT,x,T) = 



. WKMkOT.x.T)) ^ • ta »(k0T.y«b40,T«) =0.31774 

DmacCH4(x,T) 



Microscopic Thiele Modulus for C7H1 4 



jg 4n7(kOT,x,T) 



3 r 



pcatpart 



(— )- 

\Pref/ 1 



kOT 



DeffC7Hl4(x,T). 
Macroscopic Thiele Modulus for C7H14 



dpcat 
6 



DmacC7Hl4(x,T) 



♦*i7(k0T,. 5, Tin) =0.03357 



fU <W7(k0T,x,T) := If pcatpart ^/Mj.^Tl- ^^^^P^ I^ <Dmac7(k0T,ych40, Tin) =0.59265 



RatekCH4(k0T,x,T) =-ijsph(*KkOT,x,T))-iiQ)K , I , inac(kOT,x,T))- pcatpart J^J kOT 



RalekCH^kOT.O.Tin) =-1610.071 13 -sec 1 



RatekC7H14(kOT,x > T) :=-i)sph(4^(kOT,x,T)) nfpK«W7(kOT,x,T))- pcatpart- ||^JkOT 



RatekC7H14(k0T,0,Tm)=-1492.98186 sec 1 



Defines transport restricted rate constant function given catalyst effective first order rate constant - 
^ kOT, methane or cyclohexane conversion - x, and surface temperature - T 



** *****! 



r#***«***+* ******** ******************* 



Ts=Tin k0T = 0.01* sec 1 -gm" 1 Tavg =952.8725 'K 

kgpb(Tm,0.1 cm) =26.77306' sec" 1 'cm RatekCH4(kOT,0,Ts) ="1610.071 13 sec" 1 Lcat =0.002 'cm 



m+ kgpb(Tavg,0.1.cm) . J^-*WW .^mM7*K Tout = 1094.817'K 



/ kgpb(Tavg,.l-cm) \ 
\ RatekCH4(k0T,0,Ts) Lcat/ 



+ 1 



h(Tavg,.lcm) RTs 



xs=ych40 xs =0.01183 Tin =810.928'K 

Rk =6.33 m" 2 sec* 1 Ratech4 := 0.036- nf 2 - sec' 1 kOT= 10-kg' 1 'sec" 1 

O 

Set up subroutine/function to solve for the surface temperature given a rate constant 
£= TOL =0.0001 ych40 = 0.01183 x =0.01 

H Ts =1673 K Ts7 =1473 K 
Vfp Given 

ff 1 y(ych40,x)„— — 

W t^t™. -AHcomb(Ts) 



T^Tin + 

h 



Lho \ 
'100/ 



1 +. 1 

P -RatekCH4(k0T,x,Ts)Ixat 
^ Tshi(k0T) =Find(Ts) Tshi( 100 kOT) = 1010.10631 K 

nj 

Given 



y7(yOc7hl4,x) -— 

T ,_ T . A ^ ^ 6 R Ts -AHcomb7(Ts) 

Ts7=Tin-H j 

1 1 hTin," 10 



RatekC7H14(k0T,x.Ts) Lcat kgpb7 [ Tin Uo\ \ 100 



100/ 

Tshi7(kOT) = Find( Ts7) Tshi7( 100k0T) = 873.32275-K 

Rk := RatekCH4(k0T,xs,Ts) Lcat Rk7 :=- RatekC7H14(k0T,xs,Ts) Lcat 

Rk= 6.72632' sec" 1 -cm Rk7 =6.4 1773* sec" 1 'cm k0T = 10-kg~ l sec 1 

P 

xs 

^ RatecM = — kgpb(Tin,Lhc) = 15.5934-sec _1 'cm Ratech4 =0.012l5 m^-sec 1 

Rk kgpb(Tin,Lhc) kgpb7(Tin,Lhc) =4.50626* sec 1 cm 



I Ts:=Tin + R»tech4- AHcomb < Ts > .ych40 D , v , Ani<M< ^ 

^ h(Tm,Lhc) Ralech4 =0.01215m "sec 

Ratech4 Rk=6.72632-sec"' 'cm 

xs '•=■ 

Rk 

Ts = 811.4626K xs =0.1 8057' m" 3 Tin =*810.928'K 

kOT = 10-kg" 1 -sec" 1 Tshi( lO^kOT) » 1010.10631'K Tshi(0.01kg 'sec l ) = 810.96346K 

Tshi7(10akOT) = 873.32275*K Tshi7(0.0I kg 1 sec" 1 ) =810.96393*K 

xn=0,l.. 30 

kOToalo =10 5 kg'^sec 1 Tsurf =TshifkOTcab ) Tsurf7 = Tshi7(kOTcak ) 

xn ° xn V xn/ xn \ xn/ 



i . a 



ru 



1050 



1000 - 



950 ~ 



Turf 



T«rf7 




IcOTulc 



This plot shows the surface 
temperature within a honeycomb 
channel as a function of the specific 
catalytic rate constant for methane 
and methytcyclohexane combustion 
at a location within 1-mm of the 
leading edge 

Takes into account the internal pore 
diffusion (micro- and macropores) and 
boundary layer heat and mass 
transport. Upper surface temperature 
limit approaches the adiabatic 
temperature for methane and a much 
lower (Prater) temperature for C7H14 

kO unit - secM per kg washcoat 
Tunit-K 



We can also determine the actual effective first order rate constants for methane, propane, and other fuels and 
estimate the temperatures required for catalytic light-off and extinction, with roughly estimated Nu numbers. 

Other factors, such as the upstream surface temperature and gas phase radial temperature profile and thermal 
conduction of the catalyst wall, can also effect the light-off and and the critical extinction points. These calculations 
are beyond the scope of the current analysis. 



RegionO * low T PdO to 625K Region 1 - Intermediate T Region2 - Flat 909K Region3 - 
ignores inhibition by H20 and C02 625K to 909K to 1 048K decomposition Pd metal 



Apdl-ApdOcxpf-^ 2 ^^ 1 
\ 625 K 625 K 



ApdO :=24-lO 3 -20O-gm" 1 -flec~ 1 

Apd3 =Apd2expf El P d3 ) 
\ 1048. IK/ 

Interpolation formula for the four PdO+Pd rate constants 



kpd(T) := 



Apd3 2 exp 



Apd2:=Apdl«cp[^L) 
'909.1 K/ 



(Apd2 2 + Apdr 2 exp(^^) + ApdO^exp 



mi 



0.5 



O nx =0,1..41 T =500 K+nx-25 K 

„ OX 

£ CATALYTIC RATE CONSTANT NEEDED FOR HEAT BALANCE IN A 
W MONOLITH WITH FULLY DEVELOPED FLOW 



h 

MB- 



PJ 

ru 
o 

y. 



Compared with extrapolated but measured 
rate constants (mol/atmCH4/s/gcat) 



kOTcilc. 



iooaK 1000-K 



Lcat = 0.00002 -m 




T„ Taurf 
nx xo 



Compared with extrapolated and estimated 

rate constants for C7H14 (10x mol/atmCH4/s/gcat) 



kOTcak 



lil 



m 




0.6 0.8 



1.2 1.4 
1000-K 100OK 
T 'Tsmf7 



0 



fy Define functions to calculate surface temperature profiles given gas temperature and composition 



m TOL =0.1 Ts = 811.4626'K 
H* Tsx = 1500K Tsx7 = 1400K 



Given 



x = 0.01 



1 



1 



- RatekCH4(kpd( Tsx) , x, Tsx) Lcat kgpb( Tin, Lhc) 



= 7.01896-m K sz 



y(ych40,x), 



Tsx=Tin + 



0 R Tsx 



1 



1 



-AHcomb(Tsx) 
h(Tm,z) 



RatdcCH4(kpd(Tsx),x,Tsx)Lcat kgpb(Tm,z) 



Tsrf(Tin,x,z) = Find(Tsx) 



Tsch4x = Tsrf(Tra,x,0.1cm) 
Tsch4x = 823.23891 "K 



kpd(Tsch4x)=3.51311-kg 'sec 1 



Tin=810.928*K 



Z inlet = 0.1 cm 
x = 0.01 Lcat =0.002 -cm 



y(ych40,x) 0 =0.01 171 



kgpb{Tin,0.1-cm) =26.77306*sec 1 'cm 



-RatekCH4(kpd(Tsch4x),x f Tsch4x)Lcat = l.l7316*sec 1 -cm 



Given 



y7(y0c7hI4,x> P 



Tsx7=Tin + 



6 R Tsx7 -AHcomb7(Tsx7) 



\ , 1 h(Tin,z) 

- RatekC7H14( 1 0kpd( Tsx7 ) , x, Tsx7 )Lcat kgpb7( Tin, z) 

Tsrf7(Tin,x,z) =Find(Tsx7) 

Tsc7hl4x =Tsrf7(Tin,x,0.1cm) Zinlet = 0.1 cm 

Tsc7hl4x=871.51735*K Tin=810.928'K x = 0.01 Lcat = 0.002 'cm y(ych40,x) 0 =0.01171 

10 kpd(Tsc7hl4x) =44.23 122-kg" 1 -sec" 1 -RatckC7H14(10kpd(Tsdi4x),x,Tsc7hl4x)Lcat =9.2 1644* sec" 1 -cm 

kgpb7(Tin,0.1cm) = 12.3599-sec~ l 'cm 

•yy Assuming the rate constant for C7H14 combustion is ten times the rate constant for CH4 combustion, the surface 

temperature near the inlet (1-mm in) is 871 .5 K (1 100F) and near the transport limit for C7H14 and 823.2K (1050F) and 
yj surface rate controlled for CH4. If this calculation is repeated for a series of steps downstream with corrections to the 
l" 5 i gas phase concentrations and temperatures, new surface temperatures and rates can be calculated throught the length 
of the monolith. 

jS-Bj ***** ******** ********»♦****♦*♦****♦»»»»* ******♦*»*»* * »♦»****«*»*♦»*♦♦»♦**♦* »»»»^»»»»»»»»«»*»»**»*«^ 
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